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Summu*y.— The  flow  in  a round  jet  issuing  from  an  orifice  in  the  same  direction  as  a general  external  stream,  is 
investigated  theoretically  as  aji_£xtension  of  the  problem  of  a jet  issuing  into  still  air.  The  flow  in  the  upstream 
part  of  the  jet'fKegion  A of  FigTjy  in  whictLaxore  of  fluid  of  uniform  velocity  exists,  and  the  flow  in  the  downstream 
or  developed  part  of  the  jet  "{Region  B of  Fig  fy,  are  investigated  separately ; the  solutions  fitted  together  at  the 
section  at  which  the  core  disappears.  The  deviation  of  the  outside  stream  due  to  inflow  into  the  jet  is  also  considered. 
Numerical  solutions  are  derived  for  several  values  of  the  ratio  of  jet  exit  velocity  to  stream  velocity.^— 
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1.  Introduction. — The  problem  of  the  flow  in  a round  jet,  when  the  air  surrounding  the  jet 
is  at  rest,  has  been  fully  investigated1*  *•  *.  The  flow  in  and  near  a jet  with  an  external  stream 
parallel  to  the  jet  axis  is  considered  in  the  present  report,  with  special  reference  to  the  inflow 
induced  outside  the  jet.  The  fluid  is  assumed  to  be  incompressible  but  the  application  to  jets 
of  heated  or  compressible  fluids  is  briefly  considered  in  para.  7. 

• The  jet  is  assumed  to  have  a uniform  velocity  over  its  cross-section  at  the  nozzle  exit.  Mixing 
with  the  outside  stream  occurs  at  the  edge  of  the  jet,  and  the  diameter  of  the  central  core  of 
fluid  of  uniform  velocity  decreases  with  distance  downstream,  the  core  eventually  dis- 
appearing altogether.  Downstream  of  the  apex  of  the  core  the  velocity  on  the  jet  axis 
begins  to  fall  and  finally  a quasi-steady  state  is  reached  for  which  the  velocity  profiles  at  all 
sections  are  similar 
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The  presence  of  the  external  stream  makes  the  problem  more  complicated  than  for  a jet  with 
the  surrounding  air  at  rest,  and  the  method  of  solution  has  been  simplified  by  specifying  a shape 
for  the  velocity  distribution  across  the  jet  at  all  normal  sections.  For  the  first  part  of  the  jet, 
in  which  a core  of  undisturbed  fluid  exists,  the  radius  of  the  core  and  the  external  radius  of  the 
jet  are  determined  from  the  equations  of  motion,  and  after  the  core  has  disappeared  the  velocity 
on  the  axis  and  the  jet  radius  are  determined  in  the  same  way.  The  two  solutions  are  then 
joined  by  making  the  velocity  and  radius  continuous  at  the  section  containing  the  apex  < f the 
core,  and  in  this  way  a complete  picture  of  the  development  of  the  jet  is  obtained. 

The  flow  in  the  mixing  region  is  assumed  to  be  turbulent  and  the  shear  stress  is  determined 
by  application  of  the  momentum  transfer  theory’.  Reichardt's  criticisms*  of  this  theory’  do  not 
apply  to  the  present  application,  as  it  is  not  used  to  determine  the  shape  of  the  velocity 
distribution  curve  ; it  would  be  possible  to  present  the  theory  in  Reichardt’s  fohn,  but  the 
momentum  transfer  theory  has  been  retained  on  account  of  its  greater  familiarity. 

2.  Notation.— The  following  notation,  illustrated  in  Fig.  I,  will  be  used  : — 


Fig.  1 — Diagram  illustrating  Notation. 


Distance  from  jet  exit  along  axis  of  jet. 

Distance  from  jet  axis. 

Radius  of  jet  at  exit. 

Radius  of  outer  boundary  of  jet. 

Radius  of  inner  boundary  of  mixing  region,  i.e.  radius  of  core. 

Stream  Velocity. 

Jet  exit  velocity. 

uju. 

Excess  velocity  on  jet  axis  over  stream  velocity  (after  disappearance  of  core). 
Velocity  parallel  to  jet  axis  at  point  (x,  r). 

Velocity  normal  to  jet  axis  at  point  (*,  r). 

Stokes’s  stream  function  for  flow  outside  the  jet. 

Jet  momentum/2*. 

Stream  density. 

Shear  stress. 

Mixing  length. 

Mixing  length  parameter. 

Defined  in  Fig.  12  {Appendix  III). 

* This  requires  that  the  radial  inflow  velocity  into  the  mixing  region  from  the  core  shall  vanish  and  » probably  not 
exactly  true. 
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3.  hlou-  in  Region  Near  Jet  Exit. — In  this  paragraph  the  flow  in  the  part  of  the  jet  which 
contains  a core  of  fluid  of  uniform  velocity  (Region  A of  Fig.  1)  will  he  considered.  At  the 
jet  exit,  which  is  a circle  of  radius  a,  the  jet  has  a uniform  axial  velocity  Uv  It  will  be  assumed 
that  the  pressure  variations  in  and  near  the  jet  can  be  neglected  and  that  the  velocity  throughout 
the  central  core  is  constant*. 


The  momentum  equation  for  the  flow  across  a circle  of  radius  r,  whose  plane  is  normal  to  the 
jet  axis  and  at  a distance  x from  the  exit,  has  the  form  (Ref.  4,  p.  133) 

(1) 


where  u is  the  mean  axial  velocity  at  the  point  (x,  r)  and  t is  the  shear  stress  in  the  axial 
direction  at  (x,  r)  due  to  turbulent  mixing. 

The  assumed  axial  velocity  distribution  in  the  mixing  region  is  given  by  the  equation 


u — U0  -f- 


(2) 


and  is  shown  in  Fig.  2,  where  U0  is  the  external  stream  velocity  and  r0,  rv  are  respectively  the 
external  and  internal  radii  of  the  mixing  region  at  the  station  x. 


Fig.  2 — Velocity  Distribution  across  Jet  in  Region  A.  (£qn.  2). 


The  usual  formula  for  the  shear  stress,  on  the  momentum  transfer  theory  of  turbulent  flow 
(Ref.  4,  p.  207),  is 

-Kv)’  •• 


and  this  formula  will  be  adopted.  Here  l stands  for  the  mixing  length,  which  will  be  assumed 
to  be  proportional  to  the  width  of  the  mixing  region  : we  therefore  put 

l=*c{r,-rj,  (4) 


where  c is  a constant,  characteristic  of  turbulent  motion  for  the  problem  under  consideration. 


Hence 


(5) 


The  quantities  rt  and  will  be  determined  by  satisfying  equation  (1),  with  u and  r given 
by  (2)  and  (5),  at  two  positions.  The  first  position  is  at  the  outer  boundary  of  the  jet  r = r*, 
where  r = 0 and  u — so  that  equation  (1)  becomes 

JL£V«(«-tM*  = o.  ' 


At 


This  equation  is  integrated  and  the  constant  of  integration  determined  from  the  conditions  at 
the  jet  exit.  We  obtain 

j — U0)  dr  — M — {Ux  — U0), (6) 

where  a is  the  radius  of  the  jet  at  the  exit.  The  constant  M is  proportional  to  the  jet  thrust, 
which  is  equal  to  2 nM  when  the  fluid  is  taken  in  from  the  stream  ; but  if  the  fluid  is  supplied 
by  other  means  the  jet  thrust  is  increased  by  the  source  thrust  of  the  added  fluid. 

The  second  position  at  which  equation  (1)  will  be  satisfied  is  half-way  across  the  mixing 

region,  where  r = -°—,r — 1 . At  this  position,  from  equation  (2), 


u — u*  4~  ^'i  — _ !L{  ~~  Uq\ 

2 ’dr"  2 V r0  - fj  /’ 


so  that  equation  (5)  gives 


a 4 


With  these  expressions  equation  (1)  becomes 

- U>  ± ¥>.)  dr=  — + rj  (7) 

The  next  step  is  to  substitute  the  form  (2)  for  the  velocity  w into  the  momentum  equations  (6) 
and  (7).  The  integrals  occurring  can  be  evaluated  without  difficulty  and  (6)  and  (7)  reduce  to 

<*n  V + 2<ii0  rt  + aM  r0*  = a*  . . . . . , . . (8) 

[An  r,!  4 2A  „ r,  ,,  4-  A„  r,»  J = B(r,  4-  rj (9) 

where  a„  - 2 - ±),  2(»+i-  - +), 


An  — 


13  + 3a 
16 

1 - l 
16 


1 ±1 

2* 


5 + 3A  ^ 1 — 1 , 5 + 31 

A.J*  » tA  ‘ A.1 


1 + 1 5 + 31 


and  1 stands  for  the  velocity  ratio  UJUV 

The  development  of  the  jet  up  to  the  disappearance  of  the  core  is  determined  by  the 
equations  (8)  and  (9)  together  with  the  conditions  r?  = rx  — a at  the  jet  exit  x — 0.  Solutions 
have  been  derived  by  the  method  given  in  Appendix  I for  values  of  1,  the  ratio  of  the  stream 
velocity  to  the  jet  exit  velocity,  of  0,  0*125,  0*25,  0*5  and  0*75.  The  results  are  given  in 
Tables  1-5  and  in  Fig.  3,  which  also  shows  the  solution  derived  by  Kuethe4 * * * 8  for  1 — 0 by  a more 
elaborate  method  ; the  agreement  with  Kuethe's  solution  is  fairly  good  and  this  gives  confidence 
that  the  approximations  introduced  are  reasonable. 


4.  Flow  in  the  Developed  Jet. — After  the  central  core  of  uniform  velocity  has  disappeared 

due  to  the  inward  spread  of  the  turbulent  mixing  region,  the  velocity  on  the  axis  starts  to  fall 

and  further  downstream  the  velocity  distribution  across  a section  of  the  jet  settles  down  to 

a steady  shape.  No  attempt  will  be  made  to  allow  for  the  variation  of  the  shape  of  the  velocity 

distribution. 


Fic.  3 — Jet  Boundaries. 


It  will  be  assumed  that  the  axial  velocity  at  all  sections  of  the  developed  jet  is  given  by  the 
formula 

I 11  — Uq  “i  2"  + cos  ji  — )■  • ’ • • * • • * • • (10) 

I where  r0  is  the  jet  radius,  and  the  velocity  on  the  axis  is  U0  + U,  which  falls  with  increase  of 

( distance  downstream.  The  velocity  distribution  given  by  (10)  is  shown  in  Fig.  4.  In  addition, 

1 the  alternative  velocity  distribution 

« = ..  (I0a) 

also  shown  in  Fig.  4,  was  considered.  The  latter  form  probably  represents  the  correct  shape 
better  than  the  former  at  large  distances  from  the  exit,  but  the  change  was  found  to  have  little 
effect  on  the  general  development  of  the  jet,  and  the  calculations  with  it  will  not  be  discussed 
further. 


The  shear  stress  due  to  turbulent  mixing  is  given  as  before  by  (3)  but  the  mixing  length  i is 
now  given  by 

* = (11) 

which  is  obtained  from  (4)  by  putting  r,  equal  to  zero. 
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The  momentum  equation  (1)  is  still  valid  and  equation  (6)  for  the  total  momentum  still  holds 
in  the  form 

/ ° 9r  «(“  ~ U0)  dr  = M, (12) 

but  with  u now  given  by  (10).  In  addition  to  satisfying  this,  we  shall  satisfy  (1)  at  the  half- 
radius r = r0/2,  at  which,  from  (10), 


and,  from  (3)  and  (11),  1.  = — 

P 4 

Substituting  these  expressions  in  (1)  gives 


Stt  7tU 

2Tq  * 


~ j y pr  dr  - (u0  + JV  pr  « dr  - - £ P<^  * r0 (13) 


Equations  (12)  and  (13),  with  u given  by  (10),  are  sufficient  to  determine  the  variation  of  r0, 
the  jet  radius,  and  U,  the  excess  velocity  on  the  axis,  with  distance  downstream.  Evaluation 
of  the  integrals  involved  leads  to  the  equations 

p t/r0*(O- 1486  £/0  + 0-0861  U)  = M,  ..  ..  ..  ..  (i4) 

U ^ x <0'0578  U°  + 0 0476  V)  + ro  ^ (0-0914  U0  + 0-0933  U)  = - 1 -234  c* U\  (15) 
Next,  r0  is  eliminated  from  (14)  and  (15),  and  a differential  equation  obtained  for  U : — 


Jy  r (0*0914  lTn  4-  0-0933  171  - (°‘ 1486  Uo  + °' 1722  ^ (0  0578  U0  4-  0-0476C/)! 
U dx  l ' 2(0- 1486  l/0  + 0-Wr0)~ 1 J 


1 -234  c*  17*/* 


If  we  put  UJU  = x,  this  equation  becomes,  o«i  reduction. 


(0-1486  U0  4-0-0S61  U)'i*. 


dz  (0-0626  z + 0-0458)  f 0-0012  z 

dx  (0- 1486  z 4-  0-0861)1'*  [ (0- 1486  z 4-  0-0861)  (04)626  z -J-  0-0458)  J 

1 ■ 234  c*  U6 

(M/p)*'* (16) 


The  integration  of  this  equation  is  greatly  simplified  if  an  approximation  to  the  term  in 
square  brackets  on  the  left-hand  side  is  taken.  This  term  is  equal  to  unity  for  z = 0 and  tends 
to  unity  as  z tends  to  infinity,  having  a maximum  value  of  1-05  at  z — 0-7.  A mean  value 
of  1- 04  has  been  adopted  for  it  and,  with  this  simplification,  equation  (16)  becomes 


dz  (z  + 0-732)  _ 7-31  e*  UQ 

dx  (* -b  0-579)1'*  ~ (Af/p)1'*  ’ 


which  integrates  to 

{*  4-  0*579)l/ *(0-667  z + 0-990)  = 


7-31  c*  Ut 

(Af/p)»/> 


x 4-  const. 


(17) 

(18) 


& M'' « rt  «***-.-***  *'**’ 


Hie  next  step  is  to  link  this  solution  with  the  solution  with  the  central  core,  already  derived 
in  pat  a.  3.  We  note  first  that,  from  (6), 


?/?}*-  _ 1 f £?,  ITr  fM  / 1 - /-V  - 

f'o  " t/0[  2 (L’  L°^j  ~ \2?:-  ) a' 


Also,  front  para.  3 and  Fig.  3,  the  distance  from  the  jet  exit  of  the  point  of  disappearance  of 
the  core  is  given  by 

fj  _ . 

a ’ 

where  b is  a number  which  depends  upon  the  velocity  ratio  A ; at  this  point  z = 

U U j — U0 

— j — -r  • Substituting  these  values  in  (18)  fixes  the  constant  of  integration,  and  this  equation 
becomes 

(^r  4-  0-579)1;'(0-667  + 0-690)  - (Tj--v-  4-  0-579)1 ’”(t~TT  + O'690) 


Equation  (19)  detemtir.es  the  variation  of  U,  the  excess  velocity  on  the  axis,  v/ith  distance 
downstream  ; when  V is  known,  the  jet  radius  r0  is  given  by  (14).  There  is  necessarily  a smooth 
join  in  the  jet  radius  and  velocity'  at  the  junction  of  the  two  solutions  because  of  the  forms  (2) 
and  (10)  adopted  for  the  \elocity  distribution'',  which  become  the  same  when  r,  is  put  equal  to 
zero  in  (2). 

i Calculations  of  the  solutions  of  (19)  and  (14)  have  been  made  for  the  same  values’*  of  A as 

t before,  viz.,  0,  0- 125,  0-25,  and  0-5,  up  to  c-xja  equal  to  0-4.  These  are  given  in  Tables  2,  3,  4 

| and  Figs.  5a,  5b  and  6,  which  show  the  jet  boundaries,  the  velocity  on  the  axis,  and  the  rate  of 

I spread  of  the  surface  on  which  the  velocity  is  half-way  between  the  stream  velocity  and  the 

j velocity  on  the  axis  at  the  same  value  of  x.  It  is  worth  noting  here  that  the  outer  boundary 

; of  the  jet,  which  is  difficult  to  define  experimentally,  may  really  be  slightly  wider  than  that 

calculated,  because  the  actual  velocity  falls  rather  more  slowly  at  the  jet  edge  than  the  calculated 
j distribution.  For  this  reason  the  " half-velocity  ” line  (Fig.  6),  should  be  more  reliable  as  an 

■ indication  of  the  jet  width. 
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Fig.  5b — Velocity  on  Axis, 
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Fig.  6 — Half-Velocity  Lines. 


5.  Determination  of  the  Parameter  c*. — Up  to  this  point  the  mixing  length  parameter  c defined 
in  equations  (4)  and  (11)  has  been  left  undetermined.  This  has  been  done  to  keep  the  analysis 
as  general  as  possible.  But,  before  proceeding  with  the  calculation  of  the  inflow  from  outside 
the  jet,  it  is  convenient  to  select  a value.  We  shall  take  c*  = 0 0067,  this  value  being  derived 
from  experimental  data  on  the  spread  of  a jet  issuing  from  a small  orifice  in  still  air  : the  anal 


of  0*00497  obtained  from  comparison  between  theory  and  experiment  for  region  A of  a jet 
issuing  into  still  air. 


y 


Aftii  the  calculations  described  below  had  been  finished,  the  above  discrepancy  in  the 
experimental  data  was  reconsidered,  and  it  is  now  thought  possible  that  different  values  of  c* 
applv  in  regions  A and  Li.  In  practice  c-  may  rise  from  about  0 005  in  region  .1  to  O' 0067  in 
region  V>,  oxer  a transition  region  downstream  of  the  core  apex. 

6 Inflow  Outside  the  Jet.— The  flow  inside  the  jet  is  determined  by  the  above  analysis  when 
r-  has  been  specified.  The  mass  flow  across  each  section  of  the  jet  can  be  calculated,  and  is 
found  to  increase  more  rapidly  than  the  increase  provided  by  the  natural  inflow  into  the  growing 
jet  from  the  external  stream  : the  jet  thus  induces  an  inflow  towards  itself.  To  do  this,  the 
pressure  inside  the  jet  must  be  lower  than  in  the  free  stream  and  this  appears  to  be  in  contra- 
diction to  the  previous  assumption  of  uniform  pressure  ; it  has,  however,  been  shown  by  Tollmien1 
that  this  reduction  in  pressure  is  extremely  small  for  a jet  issuing  into  fluid  at  rest. 

The  next  step  is  to  calculate  the  inflow  induced  at  all  points  near  the  jet.  This  is  done  in 
two  stages  ; the  first  is  based  on  the  assumption  that  the  inflow  is  purely  radial,  and  the  second 
stage  is  the  correction  of  the  first  to  allow  for  the  interaction  of  the  flow  at  neighbouring 
sections. 

If  v is  the  outward  radial  component  of  velocity  at  the  point  (x,  r)  and  y is  Stokes’s  stream 
function  for  the  motion,  we  have 

V = | r m dr,  ^ = — r v (20) 

For  region  A we  substitute  for  « from  (2)  into  this  integral  and  obtain,  for  r > r0, 


10 


tlie  above  formulae.  This  system  of  sinks  has  a strength  per  unit  length  proportional  to  t 


which  is  given  by  (21)  and  (23)  ; numerical  values  of 
in  Fig.  7,  for  the  various  jets  considered. 


fJL  are  given  in  Tables  2-4  and 
c-rt  U1  cx 


It  will  be  seen  that  there  is  a discontinuity  in  sink  strength  at  the  point  of  disappearance 
of  the  core  ; this  is  associated  with  the  discontinuity  of  the  rate  of  growth  of  the  outer  boundary’ 
of  the  jet  there,  and  is  a defect  in  the  solution  which,  however,  is  only  significant  in  the 
immediate  neighbourhood  of  the  jet  edge  at  this  particular  section. 

It  is  next  required  to  find  the  inflow  velocities*  associated  with  the  sink  distributions  shown 
in  Fig.  7.  This  is  a problem  which  is  not  amenable  to  analytical  treatment  and  it  has  been 
solved  by  breaking  up  the  jet  axis  into  lengths  over  which  the  sink  strength  can  be  taken  to  be 
linear,  and  summing  the  contributions  of  each  linear  section  at  a large  number  of  points.  The 
details  of  this  calculation  are  given  in  Appendix  III. 

Calculations  were  made  for  cs  = 0-0067  and  for  ratios  of  the  stream  velocity  to  the  jet  exit 
velocity  * = 0* 125,  0-25  and  0-5.  The  results  are  presented  in  Figs.  8-11,  which  give  contours 
of  the  angular  deviation  (in  degrees)  of  the  outer  stream  towards  the  jet  axis  ; for  other  velocity 
ratios  the  stream  deviation  at  any  point  can  be  quickly  determined  by  interpolation. 

The  results  of  similar  inflow  calculations  made  for  x = 0*25,  c:  = 0*005  are  given  in  Fig.  11 
for  comparison  with  Fig.  9 ; the  angle  of  inflow  is  slightly  decreased  by  decrease  of  c*. 

7.  Application  to  Jets  of  Compressible  or  Heated  Fluid. — The  solution  of  the  corresponding 
problem  for  a jet  of  compressible  or  heated  fluid  is  more  complicated,  as  the  energy  equation 
would  have  to  be  satisfied  as  well  as  the  momentum  equation.  It  would  therefore  be  convenient 
if  an  “ equivalent  '*  jet  of  incompressible  fluid  could  be  determined  for  which  the  inflow  was  the 
same  as  for  a real  jet  of  heated  fluid. 

Since  the  jet  momentum  M (equation  (6)  and  (12))  is  maintained  for  ail  sections,  although 
the  density,  temperature  and  velocity  of  a hot  jet  will  vary,  it  seems  probable  that  an  equivalent 
jet  of  incompressible  fluid  would  be  such  as  to  have  the  same  value  of  the  exit  momentum  M. 
as  defined  by  equation  (6). 


* Only  [the  radial  components  an  required.  The  axial  components  are  negligible  m comparison  with  the 

stream  velocity . 


Fig.  JO — Stream  Deviation  in  Degrees.  A *■  0*5,  c’*»  0-  0067. 


Flc.  1 1—Stream  Deviation  in  Degrees.  .1  — 0*25.  c*  — 0*005. 
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Using  these  expressions  integration  of  (24)  gives 


0*5  Bx  — x,r„  — gaV^g~gn  ra^(a»/ln  11  ie2) 

10  ■ " •) 


a00£tll  fl10 


t (* . ■ ..  (25) 

! *n(gn  gio)  + ^ g2gn  ^o”(gi igoo  gic2) 


To  evaluate  the  last  integral  let 


anro4ro 


ro(gn  - «io)  + Vrtnna  - V(«n«oo  - «102) 


c f r<A 

U J r„  + Vc^  - c^2  * 


where 


11  _ r 2 

\ i ll  *“ 


(«u  - «io) 


(gll  flio)_ 


and 


_ allfl00  a102 


(gn  - «io) 


2 ‘ 


Then,  rationalising  the  denominator, 


where 


We  have 


j - f ro(ro  Vcj2  — cjr^dr^  —(ir  i / \ 

1 - * J ~ ~vfl  + cj  - Cl* Co{/'  + ^ 

J - f rJL.-° r = _ f r0  Vcf-Cfy *~dr0 

J *'•*(1  + «*)  - *1*  ’ J r#*(i  + c)  - c* 

+ **>  } [2  + r,  vT4^  - c~  ” VST+^+ J*a 

“ 2u“bj  h + 


/l_  2(1  +cj 


To  evaluate  put  e,1  — c,rB*  = c,/?8*.  so  that 


k 

« 


h =--  - 

~ 

2(1+ c,) 


RJdR9 


(1  + Ci)  V - Cj2/^ 


+ C,  - <t/  Vc, 


)] 


Hence, 
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1 


2(1  + c*) 


vV  ~ Vo*)  + 


Vl  + C,  ^Vr0Vl  -f  c,  + c, 


V - c/n2)1/2  + C,V 

v'i  + c&s  - c#*yi*~cx'  _ 


With  this  value  of  I,  x is  given  as  a function  of  r0  by  (25)  and  r,  as  a function  of  r0  by  (8). 


APPENDIX  II 

Determination  of  c1  from  Measured  Rate  of  Spread  of  a Jet  in  Still  Air 

Experimental  data  obtained  from  tests  made  at  the  R.A.E.  on  the  rate  of  spread  of  a 
jet  issuing  from  a small  orifice  in  still  air  show  that  the  jet  spreads  conically.  The  cone  on 
which  the  velocity  is  equal  to  half  the  velocity  on  the  jet  axis  was  found  to  have  a semi-angle 
of  5 degrees. 

With  the  cosine  velocity  distribution  adopted  (equation  (10))  this  gives 

r<,/2  = x tan  5°,  i.e.  r0  = 0*  175x . . . . (26) 

For  zero  stream  velocity  (14)  and  (15)  become 

0 0861  pf/%2  r=  M, 

0 0476  Ul  +0  0933  U rn^-  = - 1-234  c*l/2. 
dx  ' dx 

Elimination  of  U between  these  equations  leads  to 

£ = 27-0c«. 

and  hence,  from  (26), 

c*  = 0-0065. 

The  accuracy  of  the  experimental  data  is  not  sufficient  to  distinguish  between  an  angular  spread 
of  5 deg.  and  5-1  deg.,  and  a value  of  c*  = 0*0067  has  been  adopted  as  the  calculations  were 
a little  simpler  with  this  value. 

The  approximate  form  assumed  for  the  velocity  distribution  will  give  the  jet  radius  rather 
less  than  the  actual  value,  which  is  difficult  to  specify  exactly. 


» 


APPENDIX  III 


Determination  of  Inflow  Velocity  from  the  Axial  Sink  Distribution 


As  stated  in  para.  6,  the  inflow  velocities  corresponding  to  the  sink  distributions  of  Fig.  7 
were  determined  by  breaking  up  the  latter  into  sections,  along  each  of  which  the  sink  strength 
could  be  treated  as  linear.  We  therefore  require  formulae  to  determine  the  inflow  velocity 
for  a linear  source  distribution  of  any  length  ; this  in  turn  is  easily  determined  if  the  velocity 
fields  of  (<i)  a uniform  source  distribution  of  unit  strength  between  the  points  A ( — 1,0)  and 
B (+  1,0),  and  (6)  a source  distribution  from  A to  B of  strength  proportional  to  the  distance 
from  the  origin,  can  be  calculated. 


For  case  (a)  the  stream  function  is5 


Fig.  12.  Notation  used  in  Appendix  III. 


” " f. 


(X  — i)  d(  _ _ | . j_  , __  rt*!1/2"!1 

[r*  + ix-  f)2]1'2  '[  l J J. 


= Vr*+~(*  + l)2  ~ Vr!  + (*  - 1)J 


(*  + 1)  (*- 1) 
cx  Vr * + {*  + l)s  Vr*  + (*-!)*' 


= cos  Ba  A-  cos  0B, 

where  0A  and  are  defined  in  Fig.  12. 

For  case  (6) 


so  that 


L = \'A~  [*+ix~ 

~ J Vr*  + (x  - 1)*  + Vr*  + (*  + 1)*  j'  ^ + f * Vr*+  (*  -T)2#, 


-~= T— -j-  Vr*  + (*  + 1)*  ~ Vr*  + (*—  it* 

dx  y/r*  + (*-!)*  Vr*  + (*-f-  1)*  7-  1 vr+P 


dx  y/r*  + (x  — 1)*  Vr*  + (x 
— {R4  - COS  °a)  — iRS  — cos  $M) 


••  (28) 


The  radial  velocity  v is  equal;  to ---™  so  that  — is  the  quantity  which  we  require  to 

determine.  From  (27)  and  (28)  contours  of  equal  values  of  d^jdx  for  the  two  standard  source 
distributions  were  drawn.  With  th-j  aid  of  these  the  contribution  of  any  set  of  linear  source 
distributions  to  the  value  of  &y/dar  could  be  determined  in  a straightforward  maimer  by 
summation.  - 
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TABLE  1 
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